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SUMMARY 

t  \ . 

^An  analysis  is  made  of  the  sound  produced  when  a  field  of 
vorticity  is  cut  by  an  airfoil  in  low  Mach  number  flow.  A 
general  formula  is  given  for  the  acoustic  pressure  when  the 
airfoil  is  rigid  and  the  chord  is  acoustically  compact.  This 
expresses  the  radiation  in  terms  of  an  integral  over  the  region 
occupied  by  the  vorticity;  the  integrand  contains  factors 
describing  the  influence  of  the  thickness,  twist  and  camber  of 
the  airfoil.  Explicit  analytical  results  are  derived  for  the 
case  of  a  rectilinear  vortex,  having  small  core  diameter  and 
finite  axial  velocity  defect,  which  is  "chopped"  by  an  airfoil  of 
large  aspect  ratio.  The  acoustic  signature  generally  comprises 
two  components  associated  respectively  with  the  axial  and 
azimuthal  vorticity,  the  latter  being  determined  by  the  velocity 
defect  distribution  within  the  core.  Sound  is  generated 
predominantly  when  the  core  is  in  the  neighborhoods  of  the 
leading  and  trailing  edges.  The  contribution  from  the  trailing 
edge  is  usually  small,  however,  because  of  destructive 
interference  between  sound  produced  by  edge-diffraction  of  near 
field  energy  of  the  vortex  and  that  produced  by  vorticity  shed 
into  the  wake  of  the  airfoil  in  order  to  satisfy  the  unsteady 
Kutta  condition  that  the  pressure  and  velocity  should  be  bounded 
at  the  edge.  V7hen  the  shed  vorticity  is  assumed  to  convect  at 
the  same  mean  stream  velocity  as  the  impinging  vortex,  the 
interference  is  predicted  to  be  complete,  and  no  trailing  edge 
sound  is  generated.  If  the  shed  vorticity  is  taken  to  convect  at 
a  reduced,  "near  wake"  velocity,  which  might  be  appropriate  for 
small  scale  structures  comparable  in  size  to  the  diameter  of  the 
vortex  core,  a  small  but  non-negligible  pressure  pulse  is 
radiated  from  the  trailing  edge.  A  tentative  comparison  with 
experiment  appears  to  confirm  the  presence  of  this  trailing  edge 
pulse. 
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1.  INTRODUCTION 

An  airfoil  which  is  an  accelerated,  subsonic  notion  through 
a  nominally  quiescent  fluid  produces  sound  because  of  the  un¬ 
steady  forces  exerted  on  the  fluid.  These  forces  behave  as 
acoustic  sources  of  dipole  type,  with  dipole  axis  coinciding  with 
the  line  of  action  of  a  force.  Thus,  the  radiation  from  a 
propeller  is  associated  principally  with  a  system  of  dipoles 
whose  axes  are  normal  to  the  plane  of  rotation;  the  dipoles 
correspond  to  the  distributed  lift  (or  "thrust")  forces 
experienced  by  the  propeller  blades,  each  of  which  is  accelerat¬ 
ing  in  a  circle  (Gutin  1948).  The  lift  (and  drag)  forces  also 
vary  when  the  airfoil  encounters  a  region  of  unsteady  fluid 
motion,  either  turbulent  or  organized.  In  an  unbounded  medium, 
the  unsteadiness  is  indicative  of  the  presence  of  vorticity  in 
the  vicinity  of  the  airfoil,  and  the  impinging,  unsteady  motion 
is  usually  called  a  gust.  The  lift  exerted  on  an  airfoil  by  a 
time-harmonic  gust  in  incompressible  flow  was  first  examined 
theoretically  by  von  Karman  and  Sears  (1938). 

There  are  many  practical  problems  in  which  the  gust  can  be 
atttributed  to  a  relatively  simple  distribution  of  vorticity. 

For  example,  a  tip  vortex  shed  from  a  blade  of  the  main  rotor  of 
a  helicopter  can  impinge  on  the  following  blade,  producing  a 
characteristic  noise  often  referred  to  as  "blade  slap"  (Schmitz 
and  Yu  1986).  This  phenomenon,  in  which  the  vortex  and  the 
following  blade  are  approximately  in  parallel  planes  of  mean 
motion,  has  been  examined  theoretically  by  Martinez  and  Widnell 
(1983),  Lee  and  Roberts  (1985),  and  by  Hardin  and  Lamkin  (1987). 
Similarly,  when  the  helicopter  is  in  forward  flight,  the  main 
rotor  tip  vortices  can  be  sucked  through  the  tail  rotor  disc.  In 
that  case,  the  axes  of  the  ingested  vortices  are  almost  normal  to 
the  disc  plane,  and  sound  is  produced  by  the  multiple  "chopping" 
of  the  vortex  cores  by  the  blades  of  the  tail  rotor.  According 
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to  Leverton,  Pollard  and  VJills  (1977),  the  sound  produced  in  this 
manner  can  often  dominate  the  main  rotor  and  engine  noise  during 
helicopter  landing  maneuvers. 

Schlinker  and  Amiet  (1983)  and  Ahmadi  (1986)  have  recently 
undertaken  similar,  but  independent  experimental  studies  of  the 
sound  produced  by  the  chopping  of  a  vortex.  In  both 
investigations,  a  nominally  rectilinear  vortex  was  generated  in 
an  open-jet  wind  tunnel,  with  the  axis  of  the  vortex  parallel  to 
the  flow,  by  vorticity  shed  from  the  tip  region  of  a  stationary 
airfoil.  The  vortex  impinged  on  a  rotor  located  downstream  whose 
disc  plane  is  orientated  normal  to  the  mean  flow.  Measurements 
were  made  of  the  impulsive  sound  produced  when  the  vortex  core  is 
chopped  by  the  rotor  blades,  and  of  the  various  geometrical  and 
aerodynamic  parameters  which  are  needed  to  quantify  the  details 
of  the  interaction  of  the  vortex  and  the  blades. 

Amiet  (1986)  has  proposed  a  theory  to  account  for  the  sound 
generated  in  these  experiments  that  is  based  on  previous  studies 
of  the  gust  interaction  problem  for  an  airfoil  modelled  by  a  flat 
plate  of  infinite  span  (see,  e.g.,  Amiet  1974,  1976).  It  was 
assumed  that  the  mean  axial  velocity  in  the  core  of  the  vortex 
was  equal  to  that  in  the  ambient,  mean  wind  tunnel  flow. 

However,  the  observations  of  Schlinker  and  Amiet  (1983)  and  of 
Ahmadi  (1986)  indicate  the  existence  of  a  velocity  defect  within 
the  core,  the  peak  axial  velocity  often  being  as  much  as  10% 


smaller  than  the  mean  wind  tunnel  velocity.  Schlinker  and  Amiet 
(1983)  argued  that  a  defect  of  this  magnitude  was  too  small  to 
affect  the  level  and  characteristic  signature  of  the  acoustic 
radiation.  On  the  other  hand,  in  the  absence  of  a  velocity 


defect  the  theoretical  calculations  (Amiet  1986)  revealed  that  no 
sound  is  produced  when  the  mean  chord  of  the  airfoil  is  at  right 
angles  to  the  vortex  axis,  although  in  these  circumstances  a 
finite  axial  flow  defect  would  still  be  expected  to  generate 
sound. 
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In  this  paper  a  general  theory  is  proposed  for  the  sound 
produced  when  a  field  of  vorticity  in  low  subsonic  mean  flow  is 
cut  by  an  arbitrary,  rigid  airfoil  of  acoustically  compact  chord 
(i.e.,  when  the  wavelength  of  the  generated  sound  greatly  exceeds 
the  chord).  The  radiation  is  expressed  directly  in  terms  of  the 
vorticity  and,  although  applicable  to  flows  with  quite  general 
distributions  of  vorticity,  the  theory  is  especially  suited  to 
problems  involving  discrete  or  localized  vortex  elements.  In 
particular,  the  influence  of  an  axial  velocity  defect  within  the 
core  of  a  rectilinear  vortex,  which  determines  an  azimuthal 
component  of  vorticity  within  the  core,  can  be  treated 
theoretically  in  a  manner  entirely  analogous  to  the  treatment  of 
the  conventional  axial  vorticity  component.  The  details  of  this 
case  will  be  worked  out  for  vortex-airfoil  interactions  in  free 
space,  but  the  method  to  be  described  is  readily  extended  to 
problems  involving,  say,  the  interaction  of  helical  tip  vortices 
in  duct  flows  with  downstream  control  surfaces. 

There  are  distinct  advantages  of  confining  the  discussion  to 
airfoils  of  compact  chord.  In  the  first  instance,  there  is  a 
large  body  of  applications  in  hydroacoustics  for  which  a  general 
treatment  of  the  compressible  flow  in  the  vicinity  of  the  airfoil 
would  be  unnecessarily  complicated.  Secondly,  it  enables  the 
radiation  to  be  expressed  analytically  in  closed  form  and  the 
various  different  contributing  sources  to  be  identified.  In 
addition,  the  role  of  the  unsteady  Kutta  condition  (that  the 
velocity  and  pressure  should  remain  finite  at  the  trailing  edge 
of  the  airfoil),  and  the  significance  of  the  near  wake  convection 
velocity  can  be  examined  in  detail.  It  will  be  shown,  for 
example,  that,  for  a  flat  airfoil  of  uniform  chord  and  large 
span,  the  sound  produced  by  vorticity  shed  from  the  trailing  edge 
in  order  to  satisfy  the  Kutta  condition,  interferes  destructively 
with  the  "direct"  sound  generated  by  the  impinging  vorticity  when 
near  field  energy  of  the  latter  is  diffracted  at  the  trailing 
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edge.  When  the  shed  vorticity  is  assumed  to  convect  at  the  same 
mean  stream  velocity  as  the  impinging  vorticity,  the  destructive 
interference  is  total,  and  no  sound  emanates  from  the  trailing 
edge  region.  When  the  shed  vorticity  convects  at  a  smaller, 
near-wake  velocity,  however,  the  trailing  edge  radiation  is 
small,  but  finite. 

The  aerodynamic  sound  problem  for  a  general  airfoil  of 
compact  chord  is  formulated  and  solved  in  Section  2.  In  Section 
3  the  results  are  specialized  to  a  thin  airfoil  of  essentially 
constant  chord  and  large  aspect  ratio  at  zero  angle  of  attack, 
and  the  impinging  vorticity  is  assumed  to  convect  in  a  frozen 
pattern  at  the  velocity  of  the  mean  stream.  Analytical 
representations  of  the  acoustic  radiation  produced  by  a 
rectilinear  vortex  with  axial  velocity  defect  are  derived  in 
Sections  4  and  5.  In  Section  6  a  tentative  comparison  is  made  of 
predictions  with  the  vortex-rotor  interaction  noise  measurements 
of  Schlinker  and  Amiet  (1983). 


4 


jCmjCmjCmJLm 


Report  No.  6710 


BBN  Laboratories  Incorporated 


2.  GENERAL  REPRESENTATION  OF  SOUND  PRODUCED  BY  VORTICITY  IN  THE 
VICINITY  OF  AN  AIRFOIL  OF  COMPACT  CHORD 

2.1  Aerodynamic  Sound  in  Low  Mach  Number  Flow 

Consider  the  generation  of  sound  by  an  inhomogeneous  field 

of  vorticity  convected  in  low  Mach  number,  high  Reynolds  number 

flow  past  a  stationary,  rigid  airfoil.  The  characteristic 

wavelength  of  the  sound  is  assumed  to  be  large  relative  to  the 

chord  of  the  airfoil.  Let  the  fluid  have  mean  density  p  and 

o 

sound  speed  c,  and  let  the  mean  flow  be  at  speed  U  in  the  posi¬ 
tive  Xj-  direction  of  a  rectangular  coordinate  system  (x1,x2, 
x3).  The  origin  of  coordinates  is  located  at  some  convenient 
point  within  the  airfoil;  the  x3~axis  is  taken  in  the  spanwise 
direction,  and  the  x2~axis  is  vertically  "upwards,"  at  right 
angles  to  the  mean  flow,  as  illustrated  in  Figure  1. 

V7hen  the  Mach  number  M  =  U/c  is  sufficiently  small  relative 
to  unity  that  the  convection  and  scattering  of  sound  by  the  mean 
flow  can  be  ignored,  the  generation  of  sound  by  the  unsteady  flow 
is  governed  by  the  solution  of  the  inhomogeneous  wave  equation 

{ a2/c2 at2  -  v 2 } b  =  div(urv)  ,  (2.D 

(Howe  1975),  where  t  denotes  time,  v_  -  v_(  x_,  t )  the  velocity, 
w  =  curl  _v  the  vorticity,  and  B  is  the  stagnation  enthalpy 
defined  in  uniform,  homentropic  flow  by 

B  -  /  dp/p  +  v2/2  ,  (2.2) 

p  being  the  pressure  and  p  =  p(p)  the  local  density.  At  large 
distances  from  the  airfoil,  in  the  acoustic  far  field,  the 
perturbation  values  of  the  pressure  p  (x_,  t )  and  the  stagnation 
enthalpy  B'(_x^,t),  say,  are  related  by 
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B'(x,t)  =  p(x,t)/pQ  ,  M  <<  1  .  (2.3) 

Equation  (2.1)  is  to  be  solved  for  B  in  terms  of  the 
aeroacoustic  source  w*v  ,  whose  properties  are  specified,  or  have 
been  determined  by  prior  analysis.  The  solution  must  satisfy  the 
radiation  condition  of  outgoing  waves  at  large  distances  from  the 
airfoil,  and  comply  with  the  rigid  boundary  condition  that  the 
normal  velocity  on  the  airfoil  vn,  say,  vanishes. 


FIGURE  1.  CONVECTION  AND  SHEDDING  OF  VORTICITY  IN  MEAN  FLOW  PAST 
THE  AIRFOIL. 
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A  formal  representation  of  the  solution  can  be  written  down 
by  introducing  a  Green's  function  which  is  the 

solution  of  (2.1)  with  outgoing  wave  behavior  when  the  right  hand 
side  is  replaced  by  the  impulsive  point  source  6  ( x-y.)  6  (t-T  ) .  If 
G  is  defined  to  have  vanishing  normal  derivative  3G/3xn  on  the 
airfoil  (and  therefore,  by  reciprocity,  such  that  3G/3y  =  0, 

where  xn,  yn  are  directed  into  the  fluid)  it  follows  by  routine 
application  of  Green's  theorem  (Stratton  1941)  that 

B  ( x ,  t )  =  -  £  G  ( x ,  y ;  t ,  t  )  <  (y,i)  +  UAv)  (y,x)>  dS(y)dx 

S  -  -  ryn  -  -  -  n  -  ) 


-  /  (<jj*v)(y,x)  •  —  ( x ,  y ;  t ,  t )  d3ydx  .  (2.4) 

—  —  dy  —  —  — 

In  this  expression  the  first  integral  on  the  right  hand  side  is 
taken  over  the  surface  S  of  the  airfoil  (with  surface  element 
dS(x))/  the  normal  component  (<o*v)n  being  directed  into  the 
fluid.  The  second  integral  is  over  the  region  of  space  within 
which  u  *  (). 

To  simplify  (2.4)  we  note  that  in  high  Reynolds  number, 
homentropic  flow  the  momentum  equation  can  be  expressed  in  the 
form 


3v/3t  +  VB  +  ufv  =  0  (2.5) 

Since  vn  =  0  on  the  airfoil,  it  follows  that  the  surface  integral 
in  (2.4)  may  be  discarded,  yielding 

B(x,t)  =  -  /  (uTv)(y.,t)  .  —  ( x,  y  ;t ,  x  )d3ydx  .  (2.6) 

The  volume  integral  is  taken  over  fluid  elements  where  w  *  0. 

This  comprises  the  region  containing  the  incident  vorticity  and 
that  occupied  by  vorticity  shed  from  the  airfoil.  The  latter 
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includes  vorticity  shed  from  the  trailing  edge  in  response  to 
excitation  by  the  impinging  vorticity,  together  with  tip  vortices 
when  the  airfoil  is  at  a  finite  angle  of  attack.  In  regions 
where  u  =  0^  equations  (2.5),  (2.6)  imply  that  the  unsteady 
velocity  component  is  equal  to  where  the  potential  $ 

satisfies  3<|>/3t  =  -B. 

2.2  The  Low  Frequency  Green's  Function 

In  aerodynamic  sound  problems  involving  the  interaction  of 
low  Mach  number  vortical  flow  with  a  solid  body,  the  wavelength 
of  the  principal  ("dipole")  component  of  the  generated  sound 
greatly  exceeds  the  length  scale  of  variation  of  the  geometric 
features  of  the  body.  These  variations  may  be  regarded  as 
surface  inhomogeneities  at  which  a  small  fraction  of  the  near 
field  ("hydrodynamic")  energy  of  the  vortex  system  is  diffracted 
into  sound.  In  these  circumstances  the  Green's  function 
G  (_x_,_y_;  t ,  t  )  can  be  taken  in  its  low  frequency  form  (as  discussed 
in  detail  by  Howe  1975).  In  the  case  of  an  airfoil  of  compact 
chord,  and  when  the  observation  point  x_  is  in  the  acoustic  far 
field,  this  is  given  by 

G(x,y;t,x)  =  «{t-T-|x-y|/c}  ,  |x|  ♦  »,  H  <<  1  , 

(2.7) 

where  Y^  =  Y^(x^  is  equal  to  the  velocity  potential  of  an 
hypothetical  incompressible  flow  in  the  i-direction  which  has 
unit  speed  at  large  distances  from  the  airfoil  (i.e.,  such  that 
Yi  y ^  at  large  distances).  Observe  that,  since 
3G/3yn  ~  3Y/3yn  ,  the  normal  derivative  condition  3G/3yn  =  0  on 
the  airfoil  is  satisfied  provided  the  potential  functions  Y^ 
satisfy  this  condition. 
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Making  use  of  (2.7)  in  the  integrand  of  (2.6)  we  find  (from 
(2.3))  that  the  acoustic  pressure  fluctuations  at  large  distances 
from  the  airfoil  are  given  by 

P'*'11  '  4,c|xP  |t  1  (“-vlJy-t-lx-Yl/c)  _i(y)  d*y  . 

(2.8) 

The  proportionality  of  this  result  to  the  direction  cosine  x^/|2l| 
indicates  that  the  component  of  the  radiation  which  involves 
is  generated  by  dipole  sources  orientated  in  the  i-direction. 

The  components  Ylf  Yg  account  for  sound  produced  as  a  result  of 
the  finite  thickness,  camber,  twist  and  angle  of  attack  of  the 
airfoil.  Note  that,  since  the  airfoil  has  compact  chord,  we  can 
take  _Y  =  (0,0,y3)  in  the  retarded  tine  argument  [t]  =  t  -  |x_-_y|  /c 
of  (2.8).  This  is  an  adequate  approximation  provided  that  the 
distance  from  the  airfoil  of  vorticity  at  the  retarded  position 
(X»  It)  )  of  space-time  is  much  smaller  than  a  characteristic 
wavelength  of  the  sound  it  is  producing.  But  this  is  always  the 
case  when  the  radiation  is  generated,  as  here,  by  the  motion  at 
low  Mach  numbers  of  the  near  field  of  the  vorticity  in  the 
vicinity  of  the  airfoil. 

2.3  Modeling  the  Convection  of  Vorticity 

Equation  (2.8)  expresses  the  acoustic  pressure  in  terms  of 
the  retarded  distributions  of  vorticity  and  local  convection 
velocity  v^  In  low  Mach  number  flows  the  functional  forms  of 
these  quantities  may  be  determined  to  a  sufficient  approximation 
by  analytical  or  numerical  integration  of  the  equations  of 
incompressible  flow.  Hardin  and  Lamkin  (1987),  for  example,  have 
used  a  numerical  integration  scheme  to  treat  an  analogous  problem 
of  sound  production  by  convected  discrete  vortex  elements.  For  a 
thin  airfoil  at  small  angle  of  attack,  a  first  approximation 
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would  be  to  assume  that  vorticity  is  swept  past  the  airfoil  in  a 
frozen  pattern  at  the  velocity  of  the  undisturbed  mean  flow. 
According  to  Goldstein  and  Atassi  (1976)  this  procedure  yields 
reasonable  predictions  when  the  hydrodynamic  wavelength  of  the 
vorticity  greatly  exceeds  the  thickness  of  the  airfoil. 

When  a  rectilinear  vortex  is  chopped  by  an  airfoil  moving 
parallel  to  a  plane  whose  normal  is  nearly  aligned  with  the  axis 
of  the  vortex,  the  length  scale  of  the  vorticity  in  the  chordwise 
direction  is  equal  to  the  vortex  core  diameter.  This  is 
typically  much  smaller  than  the  chord,  and  finite  values  of  the 
camber,  twist,  thickness,  and  angle  of  attack  would  then  be 
expected  to  result  in  the  distortion  of  the  vortex,  and  perhaps 
also  in  the  generation  of  bending  waves  which  propagate  along  the 
vortex  axis  (Leibovich,  Brown  and  Patel  1986).  This  could  have  a 
significant  effect  on  the  acoustic  signature,  especially  in 
radiation  directions  parallel  to  the  mean  plane  of  the  airfoil. 
Accordingly,  an  accurate  characterization  of  the  radiated  sound 
is  possible  only  after  a  detailed  investigation  of  the  dynamics 
of  the  vortex-airfoil  interaction  has  been  performed. 

In  the  following  these  difficulties  will  be  avoided  by 
adopting  the  approach  of  thin  airfoil  theory  (Ashley  and  Landahl 
1985),  and  specifically  by  modeling  the  airfoil  as  a  flat,  rigid 
lamina  at  zero  angle  of  attack.  The  tip  vortices  which  generate 
the  mean  lift  are  then  absent  and,  provided  the  mean  velocity  is 
large  compared  to  the  perturbation  velocities  induced  by  the 
impinging  vortex,  and  in  particular  relative  to  the  upwash 
component  of  velocity  (normal  to  the  plane  of  the  airfoil),  it  is 
permissible  to  assume  that  the  vortex  will  be  convected  parallel 
to  the  Xj-axis  at  the  undisturbed  mean  flow  speed  U.  Thin 
airfoil  theory  determines  the  strength  of  the  additional 
vorticity  which  is  shed  from  the  trailing  edge  during  the 
interaction  by  application  of  the  unsteady  Kutta  condition,  and 


ft 


>x> 


>*.t, 


Report  Mo.  6710 


BBN  Laboratories  Incorporated 


assumes  that  it  is  also  convected  at  speed  U.  The  length  scale 
of  the  shed  vorticity  must  be  of  the  same  order  as  the  diameter 
of  the  vortex  core.  However,  if  this  is  much  smaller  than  the 
chord  of  the  airfoil,  it  is  arguably  more  appropriate  to  require 
the  convection  velocity  in  the  wake  just  downstream  of  the 
trailing  edge  to  be  equal  to  a  near  wake  mean  velocity  Uc  <  u. 
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3.  SOUND  PRODUCED  BY  UNIFORMLY  CONVECTING  VORTICITY  IMPINGING 
ON  A  THIN  AIRFOIL  OF  INFINITE  SPAN 

3.1  The  Directly  Generated  Sound 

The  frozen  field  of  vorticity 

“  =  £(Xj-Ut , x2 ,x3 )  ,  (3.1) 

is  assumed  to  convect  at  the  uniform  mean  steam  velocity  U  and  to 
be  cut  by  an  airfoil  which  is  idealized  as  a  thin,  rigid  lamina 
at  zero  angle  of  attack.  We  then  have 


Yi  -  Vi*  y2  -  y3  -  y3  ' 


(3.2) 


where  Y2 (y )  is  the  velocity  potential  of  incompressible  flow  in 
the  2-direction  which  has  unit  speed  at  large  distances  from  the 
airfoil.  The  identity  w'v  =  (v«V)v  -  V(v2/2),  and  the  rigid 
surface  condition  vn  =  0,  then  imply  that  only  the  y2~dipole 
makes  a  finite  contribution  to  the  radiation  integral  (2.8), 
which  now  assumes  the  form 


-p  UcosG  / 

— — — —  ^  J  <n3(y,  ft]  ) 

it c  x  3t  i  d  — 


p(x,t)  = 


°  2  °  2 

_  -  n2(y,[tn— 


|  d3z  ' 


(3.3) 


where  6  =  cos-1 (x2/| x | )  is  the  angle  between  the  radiation 
direction  _x^  and  the  normal  to  the  airfoil  (the  x2-axis). 

For  an  airfoil  of  large  aspect  ratio  9Y2/3y3  <<  3Y2/3y2 
except  in  the  immediate  neighborhoods  of  the  airfoil  tips.  The 
second  term  in  the  brace  brackets  of  the  integrand  of  (3.3)  may 
then  be  discarded.  If  the  vortex-airfoil  interaction  occurs  at  a 
sufficient  distance  from  the  tips  that  the  chord  may  be  assumed 


noESEWss •  '^rv>.-«tynv»v^ww - v  •;•  -v-, v’.w'v. v. •vuvuvwijv'vi xr, jut *mxTimwTi.wvwvini* . 
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to  be  uniform  and  to  occupy  the  interval  | x x  j  <  a  of  the  x2-plane 
(see  Figure  2),  we  may  then  take 


y2  =  Re)-t(s2-a2)1/2| 


where 


=  y.  +  iy. 


} 


(3.4) 


The  use  of  this  approximation  is  tantamount  to  assuming  the 
airfoil  to  have  infinite  span. 


5 

S 


C 


-►u 


FIGURE  2, 


CONFIGURATION  OF  THE  INTERACTION  PROBLEM  FOR  AN 
AIRFOIL  OF  INFINITE  ASPECT  RATIO. 
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Let  us  next  suppose  that  the  vorticity  £  impinges  on  the 
airfoil  in  an  acoustically  compact  neighborhood  of  x3  =  0,  so 
that  the  retarded  time  may  be  taken  to  be  [t]  =  t  -  |_x^| /e . 

Define  the  Fourier  transform  n3(K),  K  =  (KlfK2,K3)  of  fi3  relative 
to  a  coordinate  system  x/,  say,  that  is  translating  uniformly  at 
speed  U  in  the  x1~direction,  i.e.,  such  that 


A  1  —  i  V  •  v  1 

°3(i>  =  TZTjT  I  «3  -  -  d3*'  , 


x'  =  x  -  (Ut,0,0)  . 


(3.5) 


Then 


03 (Xj-Ut ,x2 , x  3 ) 


1 

U 


/  n3  (u/U,K2  ,K3  )e1^-*-”“t^da)dK2dK3  , 


(3.6) 


where  we  have  introduced  the  radian  frequency  u>  =  UK1  of  the 
constituent  components  of  the  radiation  field. 


Inserting  this  representation  into  (3.3)  (with  the  second 
term  in  the  brace  brackets  discarded),  we  obtain 

3Y„ 


pn(x,t)  = 


-p  cose 
o  3 


2c  x 


at 


/  ft. (u/u,k  ,0) 


3y; 


(y ! 'Y2  ) 


x  expfi  (ujyj/U+Kjyg-u)  [t]  )}do)dK2dyjdy2  ,  (3.7) 

wherein  the  notation  pfi(x,t)  signifies  the  acoustic  pressure 
fluctuations  generated  directly  by  the  impinging  vorticity  £, 
with  no  account  taken  of  the  coherent  radiation  from  the 
vorticity  ft  ,  say,  shed  into  the  wake.  VJhen  y.  is  defined  as  in 
(3.4),  the  remaining  integrations  with  respect  to  Vl>y2  are 
readily  evaluated  (divergent  integrals  being  treated  as 
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generalized  functions,  Lighthill  1958)  by  making  use  of  an 
integral  definition  of  the  Bessel  function  J^x)  given  in 
Abramowitz  and  Stegun  (1972,  Section  9).  This  yields 


Pfl(xrt) 


irp  aUcose  «  (on,  (u)/U,K,  ,0 )  J.  (tua/u)e 

o  3  ,  3 

c | x I  at  1  u*  +  u2k2  ) 


-iu [t] 


-dK2do> . 


(3.8) 


3.2  The  Radiation  Produced  by  the  Shed  Vorticity 


The  strength  of  the  vorticity  shed  into  the  wake  is 
determined  by  application  of  the  unsteady  Kutta  condition,  that 
the  perturbation  velocity  and  pressure  are  finite  at  the  trailing 


edge.  This  will  be  satisfied  provided  that  B  and  3B/3x2  are 


finite.  Consider  first  the  contribution  to  B  from  the  Fourier 


component 


i  n3  (fpe^-'i-®^  ,  Kx  *  oj/U  , 


of  (3.6),  and  note  from  (3.8)  that,  in  calculating  the  resulting 
acoustic  radiation,  only  components  with  K3  =  0  need  be 


considered . 


-iwt 


Let  B^,(x1,x2)e  denote  the  corresponding  component  of 
the  perturbation  stagnation  enthalpy.  This  can  be  calculated 
from  equation  (2.1)  which  assumes  the  locally  two-dimensional 


form 


{32/3x2  +  3  2 /3  x2 } B^ 1  = 


(3.9) 


=  0,  {w/U,K,,0)e1^a)Xi/U+K2X2^ 


& 


w 


*  w  " w 
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where,  in  addition,  a  further  simplification  has  been  introduced 
by  noting  that  retarded  times  are  negligible  in  the  neighborhood 
of  the  airfoil.  The  solution  of  this  equation  is  given  by  (2.6) 
with  w“v  replaced  by  (0,q,0).  The  Green'  function  satisfies 


V2G  =  -6  (x-y_)6  ( t-T  )  , 


(3.10) 


with  3G/3y2  =  0  on  | y j  J  <  a,  y2  =  ±0,  and  can  be  taken  in  the 
form  (Milne-Thonson  1968,  Section  8.71) 


G  =  Re{  In  ( )  +  ln(a2/i:’-50) }  6  (x3-y3  )6  (t-T  )  (3.11) 

where  the  asterisk  denotes  the  complex  conjugate,  and 


C  =  z  +  (z2-a2)  /  ,  z  =  Xj  +  ix2  , 

1/2 

C  *  z  +  ( z2  -a2 )  ,  z  *  y  +  iy  . 

o  o  o  o  1  2 


It  follows  from  these  results  that 


3 ' *  *  C  +  D  (x2  -  a2)  +  finite  terms  , 

n  n  n  1 


a  +  0,  x2  =  0  , 


where  Cfl,  are  independent  of  JL*  1°  particular. 


(3.12) 


(3.13) 


3  B '  ' 
0 


o)2fi3  (oj/U,K2 ,0){  J  (ua/U)  +  iJjUa/U)} 


(x2 -a2  ) 


U)2  +  U2K2 


x,  a  +  0,  x,  =  0  . 


(3.14) 
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At  a  fixed  value  of  the  radian  frequency  w,  the  aggregate 


contribution  to  3B^/3x2  near  the  trailing  edge  is  obtained  by 


integrating  over  all  values  of  K,  : 


3  B ' 

— —  =  / 


3B '  ' 

-r— —  dK 
3x2  2 


!n3  (u/U,K2  ,0){  J  Ua/u)  +  iJj  Ua/U)  }dK, 


(Xf-a2  ) 


0)2 + (j2 


Xj  +  a  +  0,  x2  =  0  . 


(3.15) 


Turn  attention  now  to  the  trailing  edge  notion  induced  by 
the  shed  vorticity  which  lies  in  the  half-plane  xx  >  a,  x2  = 

0,  and  convects  downstream  at  speed  Uc  <  U.  As  before,  only  the 
component  makes  a  finite  contribution  to  the  induced  notion, 
and  we  may  set 


O'  («,K3)  -  y  (o),K3)H(x1-a)6(x2)ei(uXl/Uc+K3X3'u,t)  ,  (3.16) 


where , 


Q  3  =  /  £J^3(o),K3  )dK3dw  , 


(3.17) 


H(x)  (=1,  0  according  as  x  ^  0)  is  the  Heaviside  unit  function, 


and  yq (w , K3 )  is  to  be  determined. 


For  the  purpose  of  computing  the  acoustic  radiation 

attention  nay  be  confined  to  Fourier  components  with  K3  =  0.  Let 

B,e-*“t  denote  the  perturbation  stagnation  enthalpy  near  the 
s 

trailing  edge  that  is  induced  by  the  component  n',(w,0).  This  is 

S  j 

determined  by  the  solution  of 
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{  32/3x2  +  32/3x|  }B^  = 


~UcYoTjr^  H(Xi_a  )  6  (x2  )exp(  iuXj/U^  ) }  . 


(3.18) 


Now,  except  in  the  wake,  B^e-*^  =  -3$^/3t,  where  <|> ^  is  the 
velocity  potential  of  motion  generated  by  the  wake  vorticity. 

The  3-component  of  the  vortex  strength  of  an  element  of  the  wake 
of  length  dx:  is  equal  to  ^dxjexpf  iu>  (Xj/l^-t ) }  .  The  evaluation 
of  the  contribution  of  each  of  these  elements  to  the  potential 
which  must  satisfy  Kelvin's  theorem  that  the  circulation  about 
the  airfoil  and  wake  is  null,  is  a  problem  in  classical 
hydrodynamics  (Mi lne-Thomson  1968,  Section  13.50),  whose  solution 
yields 


*'  s  •?—  I  Re{i(ln(?-C)  -  In  U-a2/£  )  ]  }e1U1  (yi/Uc  t)dy1  , 


(3.19) 


where  £  =  yj  +  ty^-a2 ) 1//2 ,  x,  =  z+(  z2-a2  ) 1//2,  z  =  xx+ix2  . 

When  z  lies  in  the  vicinity  of  the  trailing  edge,  the 
integral  may  be  evaluated  to  give 


C  + 
s 


D  (x? -a2 ) 1/2  +  . 


♦  a+0,  x2  =  0  , 


where  Cg,  ds  are  constants,  and  specifically, 

3B  a  (1)  y  <  ^  % 

a  - — T/2  fHo  (u>a/u c)  +  iHl  (»a/U  )} 

3X2  4(x2-a2)  ^  ° 


x;  a  +  0,  x2  =  0  , 


where  are  Hankel  functions, 

n 


(3.20) 


(3.21) 
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It  follows  from  (3.15),  (3.21)  that  at  K3  =  0  and  fixed  w , 
the  Kutta  condition  will  be  satisfied  provided  3B'/3x,  +  3B'/3x, 
remains  finite  as  x1  -►  a  +  0,  x2  =  0.  This  will  evidently  be  the 
case  if 


Y  (oj/0) 

o 


-4u{ Jg  (ua/U)  +  iJj(uja/U)}  « 

{Hq 1 } (ua/U  )  +  iH j 1 ^ (ua/U  )}  — 


£1 3  ( oj/U  ,K2 ,0)dK. 
UT-+  uTk? 


(3.22) 


The  acoustic  pressure  fluctuations  ps  (x_»  t )  produced  by  the 
wake  at  large  distances  from  the  airfoil  may  be  calculated  from 
(3.3),  (3.4)  by  replacing  03  by  ft  ^  and  the  convection  velocity  U 

(displayed  explicitly  on  the  right  hand  side  of  (3.3))  by  Uc» 

The  result  can  be  expressed  in  terms  of  yq(«»>,0)  by  making  use  of 
(3.16),  (3.17): 

irp  au  cose  »  ,  ,  .  .  i 

,  ,  °  c  3  r  ,  .  -iu[t]_ 

p  (x,t)  =  - - — r— —  j  y  (u ,0  )H  (ua/U  )e  du  , 

s  -  4c  x  3 1  o  1  c 


(3.23) 


i .e. ,  using  (3.22), 

-np  aU  cose  < 

,  t .  o  c  3  , 

P  (x,t)  =  - 1 — | -  —  / 

s  -  c  x  3 1 


u>n3  (u/U,K2,0)Hj  '(ua/U  ) 
(u2  +  U2K2 ) 


{ J0  (ua/U)  +  iJx  (oia/U)}e'lai  [t] 
{HSl)(ua/Uc)  +  iHS1)(ua/Uc)} 


dK2du  .  (3.24) 


This  is  a  generalization  of  the  formula  that  obtains  on  the 
basis  of  conventional  thin  airfoil  theory,  in  which  one  takes 
Uc  =  U.  Making  the  latter  substitution  in  (3.24),  combining  the 
resulting  expression  with  the  direct  radiation  Pfl(x,t)  given  by 


v./.-.AV.V.-.V 


y 

Ci 
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(3.8),  and  making  use  of  a  well  known  Wronskian  identity  for 
Bessel  functions  (Abramowitz  and  Stegun  1972,  Section  9.1.17),  we 
obtain 


p(Xft)  = 


Pn(x,t)  +  p  (x,t) 


2ip  U2cos0  < 

o  3  r 

c  I  x  I  3 1 


'■  —  i  (ji  [  t  ] 

n3  (oi/U,K2 ,0)e  dK0du) 

(a)2  +  U2K2  ){Hg1)  ( wa/U )  +  iHj*'  (ua/U)} 

(3.25) 


This  expression  for  the  net  acoustic  radiation  is  applicable  when 
the  wake  vorticity  is  assumed  to  convect  at  the  main  stream 
velocity  U,  and  involves  the  Sears  function 

S  ( x )  =  2/hx{Hq1)(x)  +  iH^fx)}  , 

which  arises  in  classical,  incompressible  thin  airfoil  theory 
(Sears  1940). 
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3.3  Summary  of  Results 

For  ease  of  reference  we  now  collect  together  the  formulae 
derived  above  for  the  two  principal  components  of  the  vortex- 
airfoil  interaction  radiation: 


The  direct  radiation  p^(x,t)z 


p  (*,t)  = 


irp  aUcose 
o  a 


C  X 


—  / 

i\  1.  J 


at 


(i)0  j  { w/u  /  Kj  /  0  )  J  j  (wa/Uje 
U4  +  U2K2  ) 


-ius  [t] 


dK2da)  . 


(3.26) 


The  wake  generated  sound  p  (x,t): 

s  — 


.(1  ) 


-irp  aU  cose  »  u>n,  (u/U,K,  ,0)H.  (ua/u  ) 

_  .  .  ,  _  oc  3  ,  3  2  1  c 

ps  '  c | x |  at  ^  (u2  +  u2k2) 


{ J  0  (ua/U)  +  iJjUa/UjJe 


-iw  [t  ] 


{H^Ua/U  )  +  iH.(1)(a)a/U  )} 
0  c  1  c 


dKjdiu  . 


(3.27) 
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4.  ANALYTICAL  SPECIFICATION  OF  A  RECTILINEAR  VORTEX 

The  above  formulae  will  be  applied  to  the  case  of  a 
rectilinear  vortex  that  is  convected  at  the  mean  flow  velocity  U 
and  "chopped"  by  a  stationary  airfoil.  In  the  potential  flow 
region  the  vortex  has  circulation  r  and  there  is  a  finite  axial 
velocity  (defect)  within  the  core. 

Let  the  axis  of  the  vortex  be  in  the  direction  of  the  unit 
vector  n_ whose  orientation  is  specified  by  spherical  polar  angles 
with  respect  to  the  uniformly  translating  coordinate  system 
x_'  =  2L  ~  (Ut,0,0),  introduced  in  Section  3.  The  angle  0  is 
measured  from  the  positive  xj-axis,  and  41  is  the  longitude 
measured  from  the  xj  direction,  as  illustrated  in  Figure  3,  such 
that 

n  =  (n^n^ng)  =  (sin0cos<j>,  sin0sin<j>,  cos0 )  .  (4.1) 

The  origin  O'  of  the  moving  coordinate  system  lies  on  the  axis  of 
symmetry  of  the  vortex,  translates  at  the  convection  speed  U 
along  the  x^axis  and  coincides  with  the  fixed  origin  O  at  time 
t  =  0. 

There  are  two  components  of  the  vorticity  distribution.  The 
first  is  the  axial  vorticity  £r  which  is  parallel  to  _n  and  will 
be  taken  to  be  given  by 

r  ri  exp{-(r/R)2} 

fi  =  - -5 -  (4.2) 

— r  it  r2 

where  R  is  the  nominal  radius  of  the  core  of  the  vortex,  and  r 
denotes  perpendicular  distance  from  the  uniformly  translating 
axis  of  symmetry  of  the  core.  The  Fourier  transform  nr(K)  of  £r 
is  defined  as  in  (3.5),  and  is  readily  computed  by  transforming 
from  the  integration  variables  x^ ,  x^,  xj  to  cylindrical  polar 
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FIGURE  3.  (a)  COORDINATES  DEFINING  THE  ORIENTATION  OF  THE 
IMPINGING  VORTEX. 

(b)  COMPONENTS  OF  n. 
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coordinates  whose  axis  of  symmetry  coincides  with  that  of  the 
vortex.  In  this  way  one  finds 

r  H 

nr(K)  =  -j-y  6 (n*K)exp{-(KR/2)2 }  ,  K  =  |  K  |  .  (4.3) 

Second,  in  the  presence  of  a  mean  axial  flow  within  the  core 
of  the  vortex  (corresponding  to  an  axial  velocity  defect),  there 
will  exist  an  azimuthal  component  of  vorticity  Let  V  denote 

the  maximum  axial  velocity  (defect),  and  assume  that  the  local 
axial  velocity  _v^  nay  be  taken  in  the  form 


v_A  -  Vn  exp{-(r/R)2} 


(4.4) 


so  that  the  axial  volume  flux  is  just  equal  to  ttR2V.  In  practice 
the  axial  flow  core  radius  R  in  (4.4)  will  probably  differ  from 
the  axial  vorticity  radius  of  (4.2),  but  such  refinements  will  be 
ignored. 


The  Fourier  transform  of  v^JL'  )  is  similar  in  form  to  (4.3), 


viz . , 


nVR2 

V*>  =  — 


6 ( n • K )  expf - ( KR/2 ) 2 } 


(4.5) 


from  which  it  is  immediately  deduced  that. 


£a(K)  =  iK*vA(K) 


=  6  (j2*K)exp{  -  ( KP./2  ) 2  } 


(4.6) 


For  future  reference,  we  record  the  following  particular 
components  of  (4.3),  (4.6): 
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1  "3 

0  j,  ^  ( K  ^  ,  K2  1 0  )  =  •£  5(n^Kj  +  n2K2  )exp{ -(Kj  R/2sin4i  )2  }  ;  (4.7) 


„  iVR2 ( n2 +n2 ) 

fiA3^Kl'K2'°^  =  - 4irn -  +  n2  K2  )exp{  -  ( K  j  R/2s  in$  ) 2  }  . 
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5.  VORTEX  CHOPPING  NOISE 


The  results  of  Sections  3  and  4  will  now  be  used  to  derive 
analytical  representations  of  the  acoustic  pressure  when  the 
vortex  core  diameter  2R  is  small  compared  to  the  chord  2a  of  the 
airfoil.  This  condition  is  fulfilled  by  the  experiments  of 
Schlinker  and  Amiet  (1983)  and  Ahmadi  (1986). 


5.1  Sound  Generated  by  the  Axial  Component  of  Vorticity  a 


Consider  first  the  direct  component  of  the  radiation, 
defined  by  equation  (3.26).  Substituting  for  from  (4.7), 
performing  the  integration  with  respect  to  K2  and  the 
differentiation  with  respect  to  the  time,  we  find 


prn(-'t)  = 


-ipQaUrcose  | s i n<t>  | 


4irC  |-x| . fcan"Q —  /  jj  (wa/U)exp{-(o>R/2Usin<|>  )2-ia>  [t]  }dw. 

—  —00 


(5.1) 

In  order  to  obtain  the  leading  approximation  to  this 
integral  when  R  <<  a,  observe  first  that  2R/|sin<|>|  is  precisely 
equal  to  the  minimum  interval  in  the  Xj-direction  that  contains 
the  whole  of  the  projected  area  of  the  vortex  cross-section  on 
the  x2-plane  of  the  airfoil.  In  the  absence  of  contributions 
from  the  wake  vorticity,  the  principal  radiated  pressure  fluctua¬ 
tions  are  expected  to  be  generated  when  the  projected  area  is 
swept  over  the  leading  and  trailing  edges  of  the  airfoil,  and 
will  be  characterized  by  frequencies  w  ~  U  |  s  i  n  4>  (  /R .  The  waveform 
of  the  radiated  sound  should  therefore  be  more  sharply  peaked 
when  |^|  ~  ir/2.  If  it  is  assumed  that  |sin$|  =  0(1),  which  is 
typically  the  case  when  the  core  is  "chopped,"  it  follows  that 
the  dominant  frequencies  satisfy  u>a/U  >>  1.  Introducing  the 
change  of  variable 


X  =  oiR/2U  sin$ 


(5.2) 
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in  (5.1)  we  conclude  that  the  rapid  variations  and  principal 
components  of  the  acoustic  signature  are  associated  with 
contributions  to  the  integral  from  the  region  X  *  0(1).  The 
argument  of  the  Bessel  function  becomes  2Xa/R|sin$|  >>  1  for 
X  ~  0(1).  The  integrand  may  accordingly  be  simplified  by  use  of 
the  asymptotic  approximation 

1/2 

Jj  (x)  *  sgn  (x )  {  2/ir  |  x  |  }  sin{|x|  -  ir/4}  ,  |xj  >>  1  , 

for  real  values  of  x,  following  which  one  finds 

,  .  | 3/2  .  .  . 1/2 

p  U2rcos0  sin<|>  (a/irR) 
o 

Prn (-'t )  8c | x | tan6 

where/ 


{f(a+)  -  f(a_)}  , 

(5.3) 


a±  =  (2a|sin*|/R){l  ±  U[t]/a}  , 


(5.4) 


and , 

00 

f  ( a  )  =  —  /  cosfou2  -  ir/4 }  expf -y 4  }  du  .  (5.5) 

*  0 

The  value  of  this  integral  can  be  expressed  in  terms  of  the 
modified  Bessel  functions  I±i/4  (Gradshteyn  and  Ryzhik  1980,  page 
485) : 

f  (a  )  =  |  a  |  1//2exp(-a2/8  )  { (o2/8  )  +  sgn  ( a  )  I  ^  ( a2/8  ) }  . 

(5.6) 


However,  a  more  convenient  formula,  both  for  numerical  and 
analytical  purposes,  is  obtained  by  replacing  the  Bessel 
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functions  by  their  representations  in  terns  of  the  Kummer 
confluent  hypergeometric  function  M(a,b,z),  defined  by 

M  ( a ,  b ,  z  )  =  1  +  §  TV  +  §"(■§+ 1  j  fl  + - '  (5*7) 

namely  (Abramowitz  and  Stegun  1972,  Section  9.6.47): 

(7/2)  V£»“Z 

Iv(z)  =  -fU+v)  M  ( v  +  1/2 ,  2v+l , 2z )  ,  (5.8) 

where  r(x)  denotes  the  Euler  gamma  function.  This  leads  to 
f(a)  =  exp(-a2/4)^Fn^Ty  H(l/4,l/2,a2/4)+7f^yTT  M  (  3/4 , 3/2  ,  a2/4  )j 


TT37T) 


11(1/4,1/2,-02/4)  + 


2r  ( 5/ 4 ) 


M ( 3/4 , 3/2 , -a2/4 )  . 


(  5. 9a , b ) 

The  first  of  these  expressions  is  suitable  for  numerical  work,  and 
the  second  is  obtained  by  use  of  the  Kummer  transformation 
e“zM (a , b, z )  =  M(b-a, b, z  ) . 

The  function  f(a)  >  0,  and  has  a  single  maximum  in  the 
neighborhood  of  the  origin.  Thus,  equation  (5.3)  indicates  that 
the  signature  of  the  directly  radiated  sound  prfi(x,t)  exhibits 
equal  and  opposite  pulses  respectively  at  U[t]/a  =  +1,  i.e.,  when 
the  retarded  position  of  the  core  of  the  vortex  is  in  the 
immediate  vicinities  of  the  leading  and  trailing  edges  of  the 
airfoil.  It  is  also  evident  from  (5.3)  that  the  signs  of  these 
pressure  peaks  change  as  the  polar  angle  6  made  by  the  axial 
vorticity  vector  ft  with  the  positive  x, -direction  increases 
through  ir/2  (and  that  the  radiation  is  null  when  0  =  ir/2), 
although  the  shape  of  the  pulses  depends  on  the  value  of  $  but  not 
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0.  The  double  pulse  profile  is  illustrated  by  the  dashed  curve  in 
Figure  4,  which  depicts  the  variation  of 

c  1 21 1  Prn  fc  ^ /p0u2  rcos®  » 

as  determined  by  (5.3)  for  0  =  85°,  $  =  90°,  and  R/a  =  0.1,  as  a 
function  of  the  retarded  position  U[t]/a  of  the  center  of  the 
vortex  in  the  plane  of  the  airfoil  scaled  on  the  semi-chord  a. 

These  detailed  predictions  are,  however,  unrealistic  inasmuch 
as  they  take  no  account  of  the  influence  of  induced  vortex 
shedding  from  the  trailing  edge.  The  radiation  from  the  latter  is 

A 

determined  by  equation  (3.27),  in  which  03  is  given  by  (4.7),  and 
an  explicit  representation  of  its  contribution  when  R  <<  a  can  be 
derived  in  the  manner  described  above,  leading  to 


p  (x,t) 
r  s  — 


p  UU  TCOS0  |sin+  I  3/2(a/irR)  1/2f  (a  ) 
o  c  — 

8c  x  tane 


(5.10) 


Evidently,  this  is  significantly  different  from  zero  only  when  the 
retarded  position  of  the  vortex  is  close  to  the  trailing  edge 
( i .e. ,  when  a  *  0  ) . 

V 


The  net  acoustic  pressure  is  obtained  by  combining  this 
result  with  the  representation  (5.3)  of  the  direct  radiation, 
which  yields: 


Pr(x,t)  » 


,  “I  •  I  3/2  ,  1/2 

p  U2rcos0  sin<|>  ( a  /it  R ) 

o  1  1 

8c  x  tan0 


{ f (o  )  -  (1-U  /U ) f ( a  )} 
+  c  - 


(5.11) 

This  formula  reveals  that,  when  the  wake  convection  velocity  Uc  is 
equal  to  the  main  stream  velocity  U,  the  radiation  produced  by  the 
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wake  exactly  cancels  the  directly  generated  acoustic  pulse  fron 
the  trailing  edge.  In  that  case,  the  pressure  signature  has  a 
single  peak  which  is  generated  as  the  vortex  core  sweeps  across 
the  leading  edge  of  the  airfoil,  as  illustrated  by  the  solid  curve 
in  Figure  4.  The  shape  of  this  curve  is  similar  to  the  pressure 
profiles  computed  numerically  (by  Fast  Fourier  transform)  by  Amiet 
(1986).  In  practice  it  is  likely  that  Uc  <  U;  typically  in  the 
near  wake  Uc  a  0.8U.  The  directly  generated  trailing  edge  pulse 
is  then  only  partially  destroyed  by  the  shed  vorticity.  This  will 
be  discussed  further  in  Section  6.  When  Uc  -*■  0,  the  hydrodynamic 
wavelength  of  the  wake  vorticity  ultimately  vanishes,  and  the 
prediction  (5.11)  reduces  to  that  applicable  when  the  Kutta 
condition  is  not  imposed.  These  conclusions  are  analogous  to 
results  previously  given  by  Howe  (1976)  for  two-dimensional  edge 
noise  problems. 

The  widths  of  the  acoustic  pulses  increase  both  as  R/a 
increases  and  also  when  |sin<(ij  decreases  from  unity.  The  effect 
of  this  latter  variation  is  demonstrated  in  Figure  5,  which 
depicts  profiles  of  pf(2c,t)  for  R/a  =  0.1,  Uc  =  U,  e  =  85°,  and 
*  =  30°,  60°,  90°. 

5.2  Sound  Generated  by  the  Axial  Velocity  Defect  in  the 

Core  (Q _ ) 

—A 

The  substitution  of  the  component  (4.8)  of  the  azimuthal 
vorticity  into  equations  (3.26),  (3.27)  enables  integral 
representations  to  be  derived  for  the  directly  generated  and  wake 
generated  acoustic  pressures  produced  by  the  axial  velocity  in 
the  core  of  the  vortex.  When  R/a  <<  1  the  values  of  these 
integrals  can  be  estimated  by  means  of  the  procedure  described 
above.  It  is  then  found  that  the  net  acoustic  pressure  PA(_x_,t) 
can  be  expressed  in  the  form 


M 
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_  1/2  1/2 

-p  U2Vcos0sin<f>  |  sin<j>  |  (iraR) 

p  (x,t)  =  ° 

A  — 

x  ff'(a+)  +  (1-Uc/U)f ' (a_)}  ,  (5.12) 

where  a±,  f(a)  are  defined  as  in  (5.4),  (5.5),  and  f*(a)  =  3f(a)/3a. 
This  derivative  is  readily  evaluated  from  the  representation  (5.9b) 
by  making  use  of  the  relation 

~(a,b,z)  =  ^M(I+a,l+b,z) 

(Abramowitz  and  Stegun  1972,  Section  13.4.8),  yielding 
f'(a)  =  e-X-P(~a2/4.).  {fl(3/4,3/2,a2/4)/r(5/4) 

-aM ( 1/4 , 3/2 ,a2/4)/r (3/4)  -  a2M ( 3/4 , 5/2 , a2 /4 )/4 r ( 5/4 ) }  . 

(5.12) 

It  is  clear  from  (5.12)  that  the  pressure  profile  is 
independent  of  the  angle  0  between  the  velocity  defect  v  and  the 

n 

x3-axis  (in  the  spanwise  direction).  This  is  because  the  acoustic 
dipole  source  is  equivalent  to  two  equal  and  opposite  monopole 
sources  on  the  upper  and  lower  surfaces  of  the  airfoil  whose 
strengths  are  proportional  to  the  volume  flux  in  the  x2-direction  in 
the  core  of  the  vortex,  and  (except  in  the  singular  cases  0  =  0,n) 
this  does  not  depend  on  the  orientation  of  the  vortex.  The  formula 
for  the  directly  generated  sound  p  (x,t),  which  is  applicable  when 

nW  — 

vortex  shedding  from  the  trailing  edge  is  ignored,  is  obtained  by 
setting  Uc  =  0  in  (5.12).  The  corresponding  profile  of  the  pressure 
field  has  the  symmetric  form  illustrated  by  the  dashed  curve  in 
Figure  6  for  R/a  =0.1  and  $  =  90°.  The  identical,  but  time- 
reversed  double  pulses  are  generated  during  the  passage  of  the 
vortex  over  the  leading  and  trailing  edges  of  the  airfoil.  VJhen  the 
Kutta  condition  is  imposed  with  Uc  =  u,  the  trailing  edge  pulse  is 
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PIGURE  6.  PAR  PIELD  ACOUSTIC  PRESSURE  PROFILES  GENERATED  BY  THE 
AXIAL  VELOCITY  DEFECT  vA  FOR  R/a  =0.1,  +  =  90°. 
- ,  THE  DIRECT  RADIATED  SOUND: 
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cpA(x# t)/p0U2V(a/| x| )c°s?  WHEN  THE  WAKE  CONVECTION 

VELOCITY  U  =  U. 
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annihilated  by  destructive  interference  with  the  wake  generated 
sound,  and  the  acoustic  signature  has  the  form  shown  by  the  solid 
curve  in  the  Figure.  When  0  <  Uc  <  U,  the  trailing  edge  pulse  is 
only  partially  destroyed. 

As  in  the  case  of  the  radiation  produced  by  the  axial 
component  of  vorticity,  the  width  of  the  acoustic  pressure 
pulses  is  governed  by  the  magnitude  of  R|sin<f>|/a.  This  is 
illustrated  in  Figure  7  when  R/a  =  0.1,  Uc  =  U,  and  for  different 
values  of  <f> . 


!.•  k'  K*  I 
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RETARDED  POSITION  OF  VORTEX  U[tJ/« 


FIGURE  7.  PROFILES  OF  THE  ACOUSTIC  PRESSURE  cpA (x , t )/pQU2 V(a/| x | )cos6 
FOR  R/a  =  0.1,  Uc  =  U,  AND  DIFFERENT  VALUES  OF  +  . 
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6.  APPLICATION  TO  ROTOR-VORTEX  INTERACTION  NOISE 

The  sound  produced  when  the  tip  vortex  from  an  airfoil  is 
periodically  "chopped"  by  a  rotor  has  been  examined  experi¬ 
mentally  by  Schlinker  and  Amiet  (1983)  and  by  Ahmadi  (1986).  The 
experimental  configurations  of  both  investigations  are  depicted 
schematically  in  Figure  8.  A  trailing  vortex  is  produced  in  low 
Mach  number  air  flow  in  an  open-jet  wind  tunnel  by  an  airfoil  A 
set  at  a  finite  angle  of  attack.  The  vortex  encounters  a  rotor 
at  the  downstream  station  B,  whose  plane  of  rotation  is 
perpendicular  to  the  mean  flow.  Observations  are  made  of  the 
vortex-rotor  interaction  noise  by  filtering  from  the  received 
signal  the  (typically  much  stronger)  dipole  radiation  that  is 
generated  by  the  accelerating  mean  pressure  loading  forces  on  the 
rotor  blades. 

In  order  to  apply  the  theory  of  Section  2  to  this  experiment 
it  is  strictly  necessary  to  take  account  of  the  finite  span  of 
the  rotor  blades,  and  of  the  influence  on  the  higher  frequency 
components  of  the  radiated  sound  of  their  non-compact  chords. 

The  present  discussion  will,  however,  be  based  on  the  detailed 
predictions  made  in  Section  5  for  a  flat  strip  airfoil  of 
infinite  span.  Finite  span  effects  are  probably  unimportant 
provided  that,  in  the  region  where  the  vortex  is  cut  by  the 
airfoil,  the  length  scale  of  variation  of  the  chord  in  the 
spanwise  direction  is  large  relative  to  the  diameter  of  the  core. 

Let  UQ  denote  the  mean  velocity  in  the  wind  tunnel,  and 
introduce  parameters  6,  y  such  that  the  vortex  strength  r 
(defined  as  in  equation  (4.2)),  and  the  axial  flow  velocity 
defect  reference  speed  V  (of  equation  (4.4))  are  given  by 

T  =  4*aUoB  ,  V  =  -yUq  ,  (6.1) 
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TWIN— BLADED 
ROTOR 


FIGURE  8.  SCHEMATIC  EXPERIMENTAL  ARRANGEMENT  FOR  STUDYING  VORTEX 
ROTOR  INTERACTION  NOISE. 
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where  2a  is  the  chord  of  a  rotor  blade.  For  the  experimental 
arrangement  of  Figure  8,  $  and  y  are  both  positive.  Thus,  the 
maximum  axial  velocity  defect  within  the  core  (which,  of  course, 
refers  to  the  maximum  core  velocity  relative  to  the  ambient  flow 
velocity  UQ )  is  negative.  According  to  Schlinker  and  Amiet 
(1983),  y  is  typically  of  order  0.1. 


The  formulae  derived  in  Section  5  for  the  vortex-chopping 
noise  should  be  applicable  provided  that  M  =  U/c  <<  1,  where  U  is 
now  the  chordwise  component  of  the  velocity  of  the  airfoil  (i.e., 
of  the  rotor  blade)  relative  to  the  vortex  at  the  spanwise 
location  at  which  the  vortex  is  cut  by  a  rotor.  Thus,  using  the 
definitions  (6.1)  in  equations  (5.11),  (5.12),  the  acoustic 
pressure  in  the  far  field  may  be  cast  in  the  form 


p(x,t)  =  p  (x,t)  +  pA(x,t) 


(6.2) 


where , 


Pr (x,t) 


Pou2Mocos0(a/|_xJ  ) 


_  B  |sin$  |  3/2  (ira/R) 1/2  /*,_  ,  _ 
2tane  +' 


(1-Uc/U)f  ( ct_  )  }  , 


(6.3) 


PA(x»t ) 


PoU2Mocose<a/|x| ) 


si. 


n$|sin*|1/Z(irR/a)1/2{f '  (a  +  )  +  ( 1-U  /U  )  f  '  ( a_  ) }  . 


(6.4) 
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In  these  expressions  a±,  f(a),  f'(a)  are  respectively  defined  by 
equations  (5.4),  (5.9),  (5.13),  and  HQ  =  UQ/c. 

The  contribution  pr(x,t)  of  the  radiation  produced  by  the 
axial  component  of  vorticity  vanishes  identically  when  6  =  90° 
(when  the  axis  of  the  vortex  is  normal  to  the  plane  of  the  rotor 
disc),  and  varies  rapidly  with  0  in  the  neighborhood  of  this 
angle.  In  the  experiments  6  is  typically  close  to  90°,  so  that 
reliable  predictions  can  be  made  only  when  the  orientation  of  the 
vortex  with  respect  to  a  rotor  blade  is  known  precisely.  This  is 
illustrated  in  Figure  9,  which  depicts  the  individual  waveforms 

Pa'  Pr *  Pa  +  Pr' 
for 

B  =  0.02,  y  =  0.1,  Uc/U  =  1,  <j>  =  90°  , 

and  for  e  =  85°,  90°,  95°.  The  profile  of  pA  is  independent 
of  9,  but  that  of  pr  is  inverted  when  0  changes  from  85°  to 
95°.  Accordingly,  the  net  acoustic  signature  exhibits  markedly 
different  characteristics  for  relatively  small  variations  in 
vortex  orientation. 

A  tentative  comparison  of  theory  and  experiment  is  made  in 
Figure  10.  Figure  10(a)  shows  detail  of  the  vortex-airfoil 
interaction  radiated  pressure  signature  measured  at  a  downstream 
station  ( 0  <  IT  <  90°)  by  Schlinker  and  Amiet  (1983)  at  their  Test 
Condition  B.  This  corresponds  approximately  to 

B  =  0.3R/a,  y  =  0.1,  U  =  155  msec-1  .  (6.5) 

The  rotor  blade  chord  2a  =  5.7  cms ,  so  that  the  duration  of  the 
interaction  between  the  vortex  and  the  airfoil  is  approximately 
equal  to  2a/U  *  0.3  msec. 
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Corresponding  pressure  signatures  predicted  by  (6.2)  -  (6.4) 
are  illustrated  in  Figures  10(b)/  (c),  respectively  for  Uc/u  = 
1.0,  0.8.  We  have  taken  R/a  =  0.15,  which  is  similar  in 
magnitude  to  the  value  observed  experimentally,  and  $  =  90°,  e  = 
95°.  The  physical  scale  of  the  axis  defining  the  retarded 
position  U[t]/a  of  the  vortex  axis  in  the  plane  of  the  airfoil 
has  been  adjusted  to  be  in  approximate  correspondence  with  the 
abscissa  in  Figure  10(a).  In  the  experiment  the  values  of  the 
angles  <f>  and  6  which  specify  the  orientation  of  the  vortex  were 
obtained  from  photographic  records  of  smoke  streaks  in  the  core 
of  the  vortex.  The  degree  of  uncertainty  in  this  procedure  was 
not  quantified,  however,  and  the  value  of  6  estimated  by  this 
means  by  Schlinker  and  Amiet  (1983)  was  actually  104.5°.  If  that 
value  had  been  used  in  plotting  Figures  10(b),  (c),  the  contribu¬ 
tion  pA(2L'fc)  from  the  axial  velocity  defect  would  have  been  small 
relative  to  pr()c,t),  and  the  positive  peak  at  U[t]/a  =  -1  would 
have  been  very  much  smaller.  In  Figure  10(c)  the  secondary  peak 
generated  at  the  trailing  edge  of  the  airfoil  (U[t]/a  =  1)  arises 
from  the  incomplete  cancellation  of  the  directly  radiated  sound 
by  that  produced  by  the  wake.  Its  approximate  correspondence 
with  the  trailing  edge  peak  of  Figure  10(a)  suggests  that  the 
hypothesis  that  Uc  <  u  may  be  valid.  The  widths  of  the 
theoretically  predicted  pressure  peaks  are  narrower  than  in  the 
experiment.  This  is  probably  because  of  the  assumption  that  the 
airfoil  chord  is  acoustically  compact,  which  is  not  strictly 
appropriate  when,  as  in  the  experiment,  the  Mach  number  M  =  U/c 
is  as  large  as  0.45.  According  to  the  numerical  results  of  Amiet 
(1986)  non-compactness  would  be  expected  to  broaden  the  pressure 
signature . 

The  positive  pressure  peaks  generated  at  the  leading  and 
trailing  edges  of  the  airfoil  are  absent  in  the  initial 
theoretical  predictions  of  Schlinker  and  Amiet  (1983)  and  Amiet 
(1986).  These  authors  took  no  account  of  the  axial  velocity 
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defect  in  the  vortex  core,  nor  of  the  possibility  that  Uc  <  U. 

In  order  to  obtain  predictions  that  exhibit  such  peaks,  they 
postulated  that  the  azimuthal  velocity  vr  ,  say,  induced  by  the 
axial  component  of  vorticity  £r ,  will  in  practice  decay  more 
rapidly  at  large  distance  r  from  the  axis  than  is  implied  by  the 
potential  flow  formula  Vp  =  r/2irr.  Their  hypothesis  involved  the 
assumption  that  the  distribution  of  axial  vorticity  could  be 
modelled  by  the  formula 

T  n 

«r  =  -—^3-  {  [1+(R/Rj  )2]exp{-(r/R)2}  -  ( R/R  1  ) 2  }exp{  -  ( r/Rj  )  2  }  , 

(6.6) 

instead  of  (4.2),  where  Rj  is  an  adjustable  length  scale  which  is 
taken  to  be  equal  to  the  chord  2a  of  the  airfoil.  A  consequence 
of  this  modification  of  flp  is  that  the  net  flux  of  vorticity 
through  a  plane  at  right  angles  to  the  axis  of  the  vortex  is 
identically  zero,  the  direction  of  the  vorticity  vector  changing 
sign  at  r/R  =  2.6  when  R/a  =  0.15.  This  is  unlikely  to  occur  in 
practice,  however,  and  must  cast  serious  doubt  on  the  validity  of 
predictions  based  on  (6.6). 
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7.  CONCLUSION 

There  are  two  principal  components  of  the  sound  produced 
when  the  core  of  a  rectilinear  vortex  is  "chopped"  by  an  airfoil 
of  large  aspect  ratio.  These  may  be  attributed  to  the  axial  and 
azimuthal  distributions  of  vorticity,  the  latter  being  non-zero 
when  there  exists  a  finite  axial  velocity  defect  within  the 
core.  In  practice  the  sound  generated  by  the  axial  vorticity  is 
usually  dominant,  except  when  the  vortex  is  cut  by  an  airfoil 
moving  in  a  plane  perpendicular  to  the  vortex  axis.  In  that 
case,  the  sound  produced  by  the  axial  vorticity  is  effectively 
null . 

The  predicted  profile  of  the  acoustic  pressure  is  dependent 
on  the  flow  conditions  assumed  to  obtain  at  the  trailing  edge  of 
the  airfoil.  The  impinging  vortex  induces  vortex  shedding  from 
the  edge  into  the  wake,  the  strength  of  which  is  determined  by 
application  of  the  unsteady  Kutta  condition.  Sound  generated  by 
the  wake  interferes  destructively  with  that  produced  by  the 
vortex  sweeping  past  the  trailing  edge.  If  it  is  assumed  that 
wake  vorticity  convects  downstream  at  precisely  the  mean  velocity 
of  translation  of  the  impinging  vortex,  the  interference  is 
complete,  and  for  all  practical  purposes  no  sound  emanates  from 
the  trailing  edge.  Most  of  the  radiation  is  then  associated  with 
the  vortex-airfoil  interaction  at  the  nose  of  the  airfoil.  The 
axial  component  of  vorticity  generates  a  single  (compressive  or 
rarefaction)  pulse  during  this  interaction,  whereas  the  pressure 
signature  produced  by  the  azimuthal  vorticity  consists  of  a  pulse 
followed  by  a  larger  pulse  of  opposite  sign. 

When  the  characteristic  hydrodynamic  length  scales  of  fluid 
motion  in  the  trailing  edge  region  are  small  relative  to  the 
chord  of  the  airfoil,  and  this  is  the  case  when  the  diameter  of 
the  vortex  core  is  small,  it  is  likely  that  vorticity  in  the  wake 
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will  convect  at  speeds  that  are  less  than  in  the  free  stream,  at 
least  in  the  near  wake,  where  the  most  important  interactions 
with  the  airfoil  occur.  In  that  case,  complete  cancellation  of 
the  vortex-generated  sound  at  the  trailing  edge  does  not  occur, 
and  secondary  acoustic  pulses,  of  smaller  amplitude  than  those 
from  the  leading  edge,  are  produced.  This  hypothesis  finds 
tentative  support  in  the  comparison  in  Section  6  of  theoretical 
predictions  with  measurements  of  Schlinker  and  Amiet  (1983). 

Further  developments  of  the  theory  are  required  to  include 
the  influences  of: 

1.  Finite  span  of  the  airfoil.  This  is  a  straightforward 
extension  of  the  work  of  Sections  2  and  3,  whose 
principal  difficulty  will  be  the  determination  of 
appropriate  potential  functions  Yj  of  the  low  frequency 
Green's  function. 

2.  Finite  camber,  twist  and  thickness  of  the  airfoil. 
Differences  in  the  mean  flow  velocities  on  the  upper  and 
lower  surfaces  of  the  airfoil  will  introduce  a  finite 
difference  in  the  arrival  times  at  the  trailing  edge  of 
the  upper  and  lower  portions  of  the  chopped  vortex. 

This  could  produce  interesting  and  perhaps  noisy 
dynamical  events  in  the  wake,  some  of  which  have  been 
revealed  in  the  recent  flow  visualization  studies  of 
Cary  (1987). 

3.  Finite  distortion  of  the  vortex  on  encountering  the  air¬ 
foil.  This  is  important  when  the  diameter  of  the  core 
is  comparable  to  the  length  scales  of  geometric  features 
of  the  airfoil  and  of  the  mean  flow  over  the  airfoil 
(Goldstein  and  Atassi  1976).  In  addition,  vortex 
distortion  can  generate  bending  waves  and  other 
hydrodynamic  disturbances  which  propagate  along  the  axis 
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of  the  vortex.  "Vortex  breakdown"  may  also  occur  if  the 
duration  of  the  vortex-airfoil  interaction  is 
sufficiently  long. 

4.  Non-compact  airfoil  chord.  This  is  unlikely  to  be 

important  in  underwater  applications.  Amiet  (1986)  has 
analyzed  such  effects  for  a  flat  plate  airfoil  of  large 
aspect  ratio  when  the  wake  vorticity  convects  at  the 
velocity  of  the  mean  flow. 
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